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ABSTRACT 

We have reanalyzed and compared the dipoles of the 1.2 Jy and 0.6 Jy (QDOT) 
IRAS galaxy samples. We find strong indications from both samples for (a) 
significant contributions to the gravitational field that shapes the Local Group 
motion from depths up to ~ 170 h^ 1 Mpc and (b) a large-scale coherence of the 
dipole anisotropy, indications provided mainly by the fact that the differential 
dipoles of large equal volume shells are aligned with the CMB dipole and 
exhibit significant dipole signals. The two IRAS dipoles are indistinguishable 
within 50 h^ 1 Mpc and beyond ~ 130 h^ 1 Mpc while the QDOT dipole, 
having a lower flux limit, continues growing with respect to the 1.2 Jy sample 
up to ~ 100 h^ 1 Mpc in agreement with Rowan-Robinson et al (1990). 

Keywords: cosmology: observations - galaxies: distances and redshifts - 
infrared: galaxies - large scale structure of Universe 



1 INTRODUCTION 

The peculiar velocity of the Local Group of galaxies with 
respect to the Cosmic Microwave Background (CMB) 
with u LG =622 km/sec towards (l,b) = (277°, 30°) is a 
well established fact (cf. Kogut et al. 1993). The most 
probable cause for this motion as well as for the observed 
peculiar motions of other galaxies and clusters (cf. Dekel 
1997 and references therein) is gravitational instability 
(cf. Peebles 1980). This is supported by the fact that 
the gravitational dipole (acceleration) of many different 
samples of extra-galactic mass tracers is well aligned 
with the general direction of the CMB dipole (cf. Yahil, 
Walker & Rowan-Robinson 1986; Lahav 1987; Lynden- 
Bell et al 1988; Miyaji & Boldt 1990; Rowan-Robinson 
et al 1990; Strauss et al 1992; Hudson 1993; Scaramella 
et al 1991; Plionis & Valdarnini 1991; Branchini & Plio- 
nis 1996). However, what still seems to be under discus- 
sion is horn which depths do density fluctuations con- 
tribute to the gravitational field that shapes the Local 
Group motion. The largest such depth is defined by the 
dipole convergence depth, R con v, which is that depth 
where the true gravitational acceleration converges to 
its final value. The outcome of many studies, using dif- 
ferent flux or magnitude limited galaxy samples, is that 
the apparent value of R CO nv differs from sample to sam- 
ple, in the range from 40 to 100 h" 1 Mpc, with a strong 
dependence to the sample's characteristic depth. This 



probably implies that the apparent dipole convergence 
is spurious, due to lack of adequately sampling the dis- 
tant density fluctuations. Only the optical Abell/ACO 
cluster sample is volume limited out to a large enough 
depth (~ 240 h~ x Mpc) to allow a more reliable deter- 
mination of Rconv which was found to be ~ 160 h 
Mpc (Scaramella et al 1991; Plionis & Valdarnini 1991; 
Branchini & Plionis 1996). Recently, this result has been 
confirmed using X-ray cluster samples, which are free of 
the various systematic effects from which the optical 
catalogues suffer (Plionis & Kolokotronis 1998). 

If there is a linear bias relation between the cluster, 
the galaxy and the underlying matter density fluctua- 
tions, as usually assumed (cf. Kaiser 1984), then the 
galaxy dipole should also have similarly deep contribu- 
tions. In this study we reanalyse the 1.2 Jy and the 
deeper QDOT 0.6 Jy IRAS galaxy dipoles, initially in- 
vestigated by Strauss et al (1992) and Rowan- Robinson 
et al (1990) respectively, with the aim of investigating 
whether there are any such indications. 



2 IRAS GALAXY SAMPLES AND 
SELECTION FUNCTIONS 

We use in our analysis the two available flux-limited 60- 
fim IRAS samples; one limited at Su m = 1.2 Jy (Fisher 
et al 1995) and the other at Su m = 0.6 Jy (Rowan- 
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Figure 1. The expected N(r) distribution, according to the 
IRAS galaxy luminosity function, for the 1.2 Jy, the QDOT 
0.6 Jy and the PSCz 0.6 Jy IRAS samples. 

Robinson et al 1990), which has a 1 in 6 sampling rate. 
The IRAS 1.2 Jy contains 5763 galaxies with |fc| > 5° 
while the QDOT contains 2086 galaxies with |fc| > 10°. 
Note that although the two catalogues are not totally 
independent, a cross correlation revealed only 105 com- 
mon galaxies (with 88 < 0.6° and Scz < 800 km/sec). 

To estimate the local acceleration field it is neces- 
sary to recover the true galaxy density field from the 
observed flux-limited samples. This is done by weight- 
ing each galaxy by _1 (r), where the selection function, 
4>(r), is defined as the fraction of the galaxy number 
density that is observed above the flux limit at some 
distance r. Therefore 

4>(r) = 7~~T [ *(L)dL (1) 

where L m in(r) — 4nr 2 i/Sii m is the luminosity of a 
source at distance r corresponding to the flux limit Sum, 
v = 60-/im and (n g ) is the mean galaxy number den- 
sity, given by integrating the luminosity function over 
the whole luminosity range, with L min — 7.5 x 10 7 h 2 L@ 
since lower luminosity galaxies are not represented well 
in the available samples (cf. Rowan-Robinson et al 1990; 
Fisher et al 1995), and L max = 10 13 h 2 L Q . Obvi- 
ously, 4>(r) is a decreasing function of distance because 
a smaller fraction of the luminosity function falls above 
the flux limit at greater distances. 

For the QDOT sample we used the Saunders et al 
(1990) luminosity function while for the IRAS 1.2 Jy 
we used the parameterised selection function of Yahil et 
al (1991). We have verified, however, that the two se- 
lection functions are indistinguishable from each other 
when applied to the same flux limit. In figure 1 we 
present the <j>(r) of the IRAS 1.2 and 0.6 Jy samples, 
for the l-in-6 (QDOT) as well as for the unavailable 6- 



in-6 sampling rate (PSCz). It is evident that although 
the QDOT sample is deeper, it samples the galaxy dis- 
tribution more sparsely than the 1.2 Jy sample. 



3 DIPOLE CALCULATION 

We determine the peculiar acceleration of Local Group 
galaxies by measuring moments of the IRAS galaxy dis- 
tribution. The dipole moment: D = ^ <j>~ 1 {r)r~ 2 f, is 
calculated by weighing the unit directional vector point- 
ing to the position of each galaxy, with the gravita- 
tional weight of that galaxy and summing over all avail- 
able galaxies with distances greater than 5 /i _1 Mpc 
(on smaller scales the observed galaxies do not ad- 
equately represent the true distribution; cf. Rowan- 
Robinson et al 1990). Similarly the monopole term is: 

m = J2<t>~ 1 i r ) r ~ 2 - 

We then estimate the gravitational acceleration in- 
duced on the LG from the distribution of IRAS galaxies 
by: 

(cf. Miyaji & Boldt 1990; Plionis et al 1993). Using lin- 
ear perturbation theory (cf. Peebles 1980) and equation 
(g) we can relate the Local Group peculiar velocity with 
the estimated acceleration by: 

u LG (r)=/3/V 9 (r) (3) 

where /3i = Q 6 /bi and bi is the IRAS galaxy to un- 
derlying mass bias factor. 

3.1 Treatment of the IRAS galaxy data 

Due to systematic effects and biases present in the data 
we have to perform various corrections to the raw dipole 
estimates. Firstly we need to treat the excluded, due to 
cirrus emission, galactic plane. We do so by extrapolat- 
ing to these regions the data from the rest of the unit 
sphere with the help of a spherical harmonic expansion 
of the galaxy surface density field and a sharp mask (cf. 
Yahil et al 1986; Lahav 1987). Secondly, about 4% of 
the sky is not covered by the catalogue and we apply 
to these areas a homogeneous distribution of galaxies 
having the mean weight, estimated from the rest of the 
sky. Thirdly, due to discreteness effects and the steep se- 
lection function with depth we have an additive dipole 
term, the shot-noise dipole, for which we have to cor- 
rect our raw dipole estimates. Assuming Gaussianity, 
the Cartesian components of the shot noise dipole are 
equal (a x — a y = a z ) and thus a\ D — 3a 2 1D (cf. Hud- 
son 1993). Taking the coordinate system such that one 
of the shot-noise dipole components is parallel to the 
z-axis of the true dipole and we can attempt an ap- 
proximate correction of the raw dipole according to the 
following model: 
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Figure 2. Difference between the Monte-Carlo and the 
Strauss et al (1992) 1-d shot-noise estimates, for both IRAS 
1.2 Jy and QDOT samples. 



(4) 



Note that this correction model although more severe 
than the usual D\ or = D^. am — a^ D model, it provides 
qualitatively similar dipole corrections. We choose, how- 
ever, to use this model in order to be conservative and 
to obtain a sort of lower limit to the resulting dipole, as 
far as the shot-noise correction is concerned, and thus 
via eq.(^|) an upper limit to the estimated cosmological 
P parameter (see section 4.4). 

To calculate ozd we use two methods; a Monte- 
Carlo simulation approach in which we randomise the 
angular coordinates of all galaxies while keeping their 
distance, and thus their selection function, unchanged 
while the second method is the analytic estimation of 
Strauss et al (1992): a\ D ~ J2 9^ Wi"^ 1 )- Figure 2 
shows the difference (in velocity units) between the two 
shot-noise estimates. It is evident that the two methods 
give equivalent results although due to the statistical 
nature of the first method we believe that it performs 
better on large depths, where the number density of 
IRAS galaxies is very low. 

3.2 z to 3d frame correction 

The final but essential correction is to transform red- 
shifts to 3-d distances in order to minimise the so called 
'Kaiser' effect (Kaiser 1987). This effect can be under- 
stood by noting that the distribution of galaxies in red- 
shift space differs from that in real comoving space by 
a non-linear term: 



cz = H r + (v(r) - v(0)) ■ 



(5) 



where v(0)(= ulg) is the peculiar velocity of the Lo- 
cal Group and v(r) the peculiar velocity of a galaxy 
at position r. If v(r) had random orientation, then 



1000 




Figure 3. Comparison of local 1.2 Jy IRAS galaxy peculiar 
velocities provided by our model and by the full iterative 
algorithm within cz ~ 3600 km/sec (see text). The diagonal 
line corresponds to v reC ons • ? = Vbulk • f . 



J v(r) ■ f d 3 r ~ and the last term of eq.(P) is domi- 
nated by the LG term; we thus obtain that in the LG 
frame (ie., when using cz — H a r), structures in the di- 
rection of our motion appear at a redshift smaller than 
their true distance in the CMB frame and thus they will 
artificially enhance the amplitude of the gravitational 
dipole. 

However, many studies indicate that local galax- 
ies have peculiar velocities not randomly oriented but 
rather participating in a coherent flow (bulk motion) to- 
gether with the Local Group (ie., v(r) ~ v(0)) within at 
least a volume of radius ~ 5000 km/sec (cf. Lynden-Bell 
et al 1988, Dekel 1994; 1997, Strauss & Willick 1995). If 
so, it would be reasonable to evaluate the IRAS dipole 
in the LG frame, since in this case cz m H r. However, 
this is not absolutely true since there should exist also 
a velocity component due to the local, non-linear, dy- 
namics acting between nearby galaxies and/or clusters 
of galaxies. We can therefore view the galaxy peculiar 
velocities as consisting of two vector components; a bulk 
flow and a local non-linear term: 



v(r) = V bulk (r) + v n ;(r) 



(6) 
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Inserting eq.(|6|) in eq.(|5|) and assuming that v(r) • f w 
^buik(r) ■ f , ie., that the dominant component is that 
of the bulk flow Q we can use the observed bulk flow 
profile, as a function of distance, given by Dekel (1994; 
1997) and combined with that of Branchini, Plionis & 
Sciama (1996) to correct the galaxy redshifts. The zero- 
point, Vbuik(O), and the direction of the bulk flow is 
estimated applying eq.(^) at r = and assuming, due 
to the "coldness" of the local velocity field (cf. Peebles 
1988), that v ni (0) ~ Vj n / (where v in f is the LG infall 
velocity to the Virgo Supercluster). Using the average 
value from the literature, ie. v in f ~ 170 km/sec, we ob- 
tain V bu ik(0) ^ 500 km/sec towards (l,b) ~ (276°, 15°). 

We test our model by comparing peculiar velocities 
that it provides with those resulting from the full dy- 
namical algorithm (kindly provided by Dr. Enzo Bran- 
chini) which estimates, using linear theory, the gravita- 
tional acceleration at the position of each galaxy and 
then recovers the real-space galaxy distances by solving 
iteratively the generalised Hubble law of eq| (cf. Yahil 
et al 1991; Strauss et al 1992). In figure 3 we present 
this comparison for relatively local galaxies in regions 
of Sp/p < 1 (since at dense regions the non-linear com- 
ponent that we neglect in our model will dominate the 
galaxy peculiar velocity). We find a good correlation 
within cz ~ 4000 km/sec which is in fact the region 
where such corrections can affect the dipole. The cor- 
relation, at larger distances, progressively fades away 
since the bulk flow amplitude is low at such distances 
and the galaxy peculiar velocities are dominated by the 
distant local dynamics. In any case at such distances 
we do have Jv(r) ■ f d 3 r « and thus redshift space 
distortions are dominated by the LG term in eq.(^) for 
which we do indeed correct the galaxy redshifts. Note 
that we have verified that the amount of scatter seen in 
figure 3 is well reproduced from our model if we include 
a randomly oriented non-linear velocity component with 
<(v ni ■ f) 2 } 1/2 w 320 km/sec. 

We have further tested the robustness of the recov- 
ered real-space distribution by performing 200 Monte- 
Carlo simulations in which we vary Vinf (and there- 
fore also the amplitude and slightly the direction of 
Vb u ik(0)) as well as the amplitude of Vb u lk(f) for all 
r's, by randomly sampling a Gaussian having as mean 
(p) the nominal velocity values and a = 2p/3. Further- 
more, we investigate how our results change when using 
the bulk-flow direction of Lauer & Postman (1994); ie., 
(l,b) ~ (343°, 52°) with \V bu ik{r)\ = 650 km/sec for 
r < 130 h' 1 Mpc. 

Finally, we would like to point out that it so hap- 
pens that the IRAS dipole, estimated in either the LG or 
the CMB frames, which should provide a sort of upper 

* in a sense we assume that the vector average of v„;(r) ■ 
r over a whole sky distribution of galaxies is ~ 0; not an 
unreasonable assumption in the limit of dense sampling. 



and lower dipole bounds respectively, differs very little 
and thus the z to 3d frame correction does not have a 
major consequence in our main dipole results. We do 
however investigate, in section 4.2, possible systematic 
effects that could be introduced by the frame transfor- 
mation procedure in our IRAS-CMB dipole alignment 
results. 



4 MAIN RESULTS 

In figure 4a we present the two IRAS dipoles in redshift 
space. We observe that they are consistent, although 
up to ~ 50 h^ 1 Mpc the 1.2 Jy dipole is systemati- 
cally higher than the QDOT dipole. However, once we 
correct for redshift space distortions (figure 4b) the two 
corrected dipoles almost coincide within 50 h -1 Mpc, as 
they should since both samples, due to their relative low 
flux limits, are good tracers of the Local Universe. Up 
to ~ 130 ft -1 Mpc the amplitudes of the two real-space 
dipoles deviate with the QDOT being larger than the 
IRAS 1.2Jy, which is to be expected since the QDOT 
sample has a lower flux limit and thus it can 'see' the 
distant matter fluctuations better than the 1.2 Jy sam- 
ple. Beyond ~ 130 hT 1 Mpc, however, the two IRAS 
dipoles coincide again. As can be seen from figure 4, 
redshift space distortions enhance significantly the real- 
space dipole (by ~ 12%) only within ~ 50 Mpc. 
The uncertainties of the velocity field model, probed by 
the Monte-Carlo simulations discussed previously, intro- 
duce a small scatter in the real-space dipole as indicated 
by the errorbars in figure 3b. Using the Lauer & Post- 
man (1994) bulk-flow leaves unaltered our main dipole 
results with only a small (~ 40 km/sec) amplitude de- 
crease. 

4.1 Evidence for > 100 h' 1 Mpc dipole 
contributions 

Between ~ 150 — 180 h~ x Mpc there is an apparent 
amplitude bump, seen in both redshift and real-space 
IRAS dipoles. This bump is accompanied by a ~ 5° de- 
crease of the misalignment angle between the two IRAS 
dipoles and that of the CMB, the overall misalignment 
angles at r = 200 h' 1 Mpc being ~ 23° and ~ 35° 
for the 1.2 Jy and QDOT samples, respectively. These 
facts suggest that this dipole amplitude bump is not due 
to shot-noise uncertainties but rather it is an intrinsic 
effect, indicating the existence of contributions to the 
Local Group motion from such large scales. Such con- 
tributions cannot be accurately determined, however, 
from the present flux-limited samples and deeper sam- 
ples are required for such a task (see Kolokotronis et al 
1996). To further investigate these probable deep IRAS 
dipole contributions we estimate the differential dipole 
in large equal volume shells (see Plionis, Coles & Cate- 
lan 1993 for an earlier attempt in z-space). We investi- 
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Figure 4. (a) IRAS 1.2 Jy (continuous line) and QDOT 
(dashed line) dipoles in redshift-space. (b) the correspond- 
ing dipoles in real-space. The errorbars (shown only for the 
QDOT dipole) are estimated from Monte-Carlo simulations 
of the velocity field model (see text). 



Table 1. Differential 1.2 Jy IRAS dipole directions, the cor- 
responding misalignment angles with respect to the CMB 
dipole, the dipole signal to noise ratio and probabilities of 
alignment within 88 (see text for complete definition). The 
errorbars are estimated from Monte-Carlo simulations of the 
velocity field model, with the last two shells having no cor- 
responding uncertainties because we take Vb u i^(r) = for 
r > 170 h~ l Mpc (see text). 



h^Mpc N gal 


S/N 


1° 


6° 


89 C mb 


Pmc 


Pf 


5-110 


4329 


4.7±0.2 


246.0 


35.5 


26±3 


0.003 


0.051 


110-139 


390 


0.6±0.2 


311.8 


84.1 


55±17 


0.522 


0.213 


139-159 


220 


1.7±0.2 


268.9 


8.5 


22±4 


0.026 


0.036 


159-175 


128 


0.6±0.2 


307.1 


31.5 


26±7 


0.073 


0.051 


175-188 


72 


0.0 


96.6 


61.4 


89 


0.655 


0.492 


188-200 


59 


1.6 


282.8 


43.2 


14 


0.029 


0.016 



gate shell sizes ranging from 2.5 x 10 6 to 8.2 x 10 6 h~ 3 
Mpc 3 which give qualitatively similar results. In tables 
1 & 2 we present the differential dipole directions and 
misalignment angles with respect to the CMB dipole as 
well as a measure of the significance of the dipole of each 
individual shell, given by: 



Table 2. 


QDOT dipole 


results (as in 


Table 1). 




ft -1 Mpc 




S/N 


1° 


b° 


8@cmb Pmc 


Pf 


5-110 


1267 


2.3±0.2 


232.5 


30.3 


38±3 0.016 


0.011 


110-139 


210 


-0.1±0.2 


98.8 


36.6 


113±26 0.89 


0.71 


139-159 


96 


1.2±0.2 


271.1 


-12.7 


43±9 0.116 


0.134 


159-175 


66 


1.2±0.1 


272.1 


31.6 


5±2 0.004 


0.002 


175-188 


53 


-0.9 


164.4 


-73.7 


125 0.55 


0.78 


188-200 


31 


0.7 


318.4 


-15.6 


60 0.55 


0.25 



S D . 

_ raw CQS f§Q cmb ^ _ f7) 
IV a-sD 

for the case of 8V ~ 5.5 x 10 e h' 3 Mpc 3 (six shells). 
We observe that for the QDOT sample there are 3 
shells with relatively small misalignment angles and 
dipole signal to noise ratios > 1, the deepest shell being 
[159-175] h' 1 Mpc, in which 86 cmb ~ 5°. For the 1.2 Jy 
sample, which although shallower has a better sampling, 
we have small misalignment angles (S6 cm b^ 27°) in the 
same shells but also in a deeper shell (188 — 200 
Mpc). Out of these four aligned shells there are signifi- 
cant dipole contributions (S/N > 1) only in three while 
the probability that these alignments are random is ex- 
tremely low. The formal probability that two vectors 
are aligned within 88 is given by the ratio of the solid 
angle which corresponds to 89, to the solid angle of the 
whole sphere, ie., Pf(S6) = sin 2 (<5#/2). 

We can now estimate the joint probability of align- 
ment, within the observed 89 cmb , of N independent vec- 
tors, which is given by: 



(8) 



Between three shells (first, fourth and sixth) the IRAS 
galaxy correlation function is zero, due to the large dis- 
tances involved, and consequently the shells can be con- 
sidered independent. Due, however, to the vicinity and 
therefore the possible correlation between the third and 
fourth QDOT shell we will consider their joint probabil- 
ities as limits. Therefore, we have that the joint proba- 
bility of alignment between the CMB and the differen- 
tial IRAS equal-volume dipole directions (for those with 
significant dipole signal S/N > 1) is: 

-■*< 



2 x 10" 



P^ OT < 8 x 10" 



Pl 23y ~ 2 x 10" 4 . 
for the QDOT and IRAS 1.2 Jy samples, respectively. 



4.2 Test for systematic alignment errors 

The observed differential dipole alignments could in 
principle be due to errors in the correction used to re- 
cover the 3-d frame in which we measure the dipole. 
For example, if redshift errors especially at large dis- 
tances, were the sampling rate is low, were comparable 
to a significant fraction of the LG velocity, then using 
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the 3-d galaxy distances estimated from eqs (g) , (g) and 
v(0) = 622 km/sec, we could artificially produce a false 
alignment of the distant shells differential dipole with 
that of the CMB. We therefore address the question 
of which are the dipole alignments induced totally due 
to our frame transformation procedure ie., for the ex- 
treme case where there is no intrinsic dipole and no 
LG peculiar velocity. We run 5000 Monte-Carlo simu- 
lations in which we destroy the intrinsic IRAS galaxy 
dipole as well as redshift space distortions by randomis- 
ing the angular coordinates of the galaxies while keep- 
ing their distances and therefore their selection function 
unchanged. On this intrinsically random galaxy distri- 
bution we apply our z to 3d space transformation and 
then measure the artificially induced differential dipole 
alignments due to the frame transformation itself. The 
coupling between the space distortion and the selection 
function results in a non-trivial alignment behaviour. 
In the first shell we observe anti-alignments while at 
more distant shells the artificial alignment effect does 
appear. For example at the last shell the median 59 is 
~ 67° instead of 90° . However, it is impossible to create 
the observed IRAS dipole alignments if there is no true 
signal present. We quantify this by measuring the prob- 
ability, p m c(59), of observing in our Monte-Carlo simu- 
lations dipole alignments as large as the observed IRAS 
differential dipole alignments. For the shells of interest 
we find that this probability is low and comparable to 
the expected pj, which implies that the frame transfor- 
mation uncertainties cannot induce the observed IRAS 
dipole alignments (see the corresponding values of p mc 
and pf in the tables). 

We conclude that the differential dipole directions 
are not randomly oriented with respect to the CMB and 
therefore we do not only have indications for significant 
dipole contributions from large depths but also for a 
coherent anisotropy extending to these large scales. 



4.3 Possible cause of the large-scale IRAS 
dipole contributions 

It is interesting that stronge evidence exist for deep 
dipole contributions from the available galaxy cluster 
data. Contributions up to ~ 20% — 30% of the total 
cumulative optical and X-ray cluster dipole, were found 
from ~ 140 — 160 h~ l Mpc depths (Scaramella et al 
1991; Plionis & Valdarnini 1991; Branchini & Plionis 
1996; Plionis & Kolokotronis 1998). Similar coherence 
of the differential dipoles in equal volume shells was 
also found in galaxy cluster case (Plionis & Valdarnini 
1991; Plionis et al 1993). These studies have shown that 
the cause of the deep dipole contributions should be at- 
tributed mostly to the Shapley concentration, a huge 
mass overdensity located at ~ 140 h~ l Mpc in the 
general direction of the Hydra- Centaurus supercluster 




Figure 5. Contour plot of the smoothed 1.2 Jy IRAS galaxy 
distribution in real space and on the supergalactic plane. 

(Shapley 1930; Scaramella et al 1989; Raychaudhury 
1989). 

To investigate further the possible cause of the 
present IRAS dipole results we have smoothed the IRAS 
1.2 Jy galaxy distribution in a 40 3 cube with a cell size 
of 10 h" 1 Mpc using a Gaussian with smoothing radius 
equal to one cell and weighting each galaxy by (j>~ 1 . 
Due to the coupling between the selection function and 
the constant radius smoothing, we correct the resulting 
smoothed distribution for a distance dependent effect, 
which we quantified using N-body simulations (details 
will be presented in a forthcoming paper). In figure 5 we 
present the smoothed IRAS 1.2 Jy galaxy distribution 
on the supergalactic plane (of 10 ft -1 Mpc width) within 
170 Mpc. The contour step is 0.4 in overdensity 
while the 5p/p = level appears as a thick continuous 
line. 

Well known structures appear in this plot; the 
largest and most evident is the Shapley concentration 
located at (X Bup , Y sup ) ~ (—120, 60), the Perseus-Pisces 
supercluster at (X aup ,Y sup ) w (60,-40), the Coma su- 
percluster at (X aup ,Y aup ) w (—20,70), the Ursa-Major 
supercluster at (X sup , Y sup ) « (100,100), the Pisces- 
Cetus supercluster at (X sup ,Y sup ) w (50,-140) while 
the Great Attractor (Hydra-Centaurus complex ?), at 
(Xsup, Y Bup ) ~ (—30,30) appears in the foreground of 
the Shapley concentration. Furthermore, figure 5 is very 
similar to the corresponding Abell/ACO cluster density 
field (cf. Tully et al 1992; Branchini & Plionis 1996, their 
figure 3), and we therefore obtain a consistent picture, 
from both IRAS galaxy and Abell/ACO cluster data, in 
which the Shapley concentration is the most probable 
cause of the deep dipole contributions while the general 
alignment of the Great Attractor, Perseus-Pisces and 
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Shapley superclusters is most probably the cause of the 
apparent coherence of the IRAS galaxy dipole. 

4.4 Q°' 6 /bi from the IRAS dipoles 

Using the real-space dipole results and equation (0), we 
can estimate the density parameter Piras- However, 
the value obtained should be considered rather as an 
upper limit since the deep contributions to the dipole, 
for which we do have strong indications, are probably 
not fully revealed by the present samples (see Kolokotro- 
nis et al 1996). Taking into account the scatter among 
the two IRAS samples, the amplitude variations at large 
depths and the uncertainties of the velocity model used 
to recover the real-space galaxy distances, we find: 

Piras^ 0.78(±0.1) , (9) 

in agreement with the QDOT analysis of Rowan- 
Robinson et al (1990) but slightly larger, although 
within la, than the 1.2 Jy results of Strauss et al 
(1992). This value of pi implies either that Q. < 0.66 
for b IRA s — 1 or Q a ~ 1 for b IRA s~ 1-28. 

5 CONCLUSIONS 

Using a consistent analysis procedure we find that 
within 50 h" 1 Mpc, both the 1.2Jy and 0.6Jy (QDOT) 
IRAS samples, give identical dipole results while beyond 
this depth the QDOT dipole increases substantially up 
to 100 /i -1 Mpc, in agreement with Rowan- Robinson 
et al (1990). Furthermore there are significant indica- 
tions for (a) dipole contributions from depths ~ 170 h~ l 
Mpc, in agreement with other recent large-scale studies 
(cf. Plionis & Kolokotronis 1998) and (b) a coherence of 
the dipole anisotropy extending to similar depths. The 
most probable cause of these deep dipole contributions 
is the Shapley mass concentration, while of the dipole 
coherence is the general alignment, on the supergalac- 
tic plane, of the Perseus-Pisces supercluster, the Great 
Attractor and the Shapley concentration, which span a 
range of ~ 200 h' 1 Mpc. 

A similar study of the complete (6 in 6) IRAS 
0.6 Jy sample (PSCz) should give better indications of 
these results, although the overall amplitude of the ef- 
fect could be probably estimated by a deeper catalogue 
(limited at a lower flux limit). 
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